Effective medium theory of binary thermoelectrics. 
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The transport coefficients of disordered media are analyzed using direct numerical simulation and 
effective medium theory. The results indicate a range of materials parameters for which disorder 
leads to an enhanced power factor. This increase in power factor is generally not accompanied by an 
increase in the figure of merit ZT, however. It is also found that the effective electrical conductivity 
and electronic contribution to the thermal conductivity are not generally proportional to each other 
in the presence of disorder. 



I. INTRODUCTION 

There has been considerable recent interest in utiliz- 
ing nanostructure and inhomogeneity to enhance thermo- 
electric performance. Nanostructure has been predicted 
to reduce phonon thermal conductivity without reducing 
electrical conductivity 111, enhance the power factor via 
quantum confinement [2[ , and provide strongly energy- 
dependent electron scattering [J-Q - beneficial features 
for thermoelectrics. One type of nanostructure consists 
of the inclusion of metallic structures in a thermoelectric 
(i.e semiconductor) framework. These types of nanos- 
tructured materials can exhibit two unusual properties: 
unexpectedly high Seeback coefficient @ , and anomalous 
magnetoresistance Q. The anomalous magnetoresistive 
properties can be understood as a consequence of disor- 
dered transport in a random resistor network model Q. 

Motivated by these models of magnetoresistance,I re- 
visit similar models of thermoelectric properties [9| to 
see if the same mechanism can explain both phenom- 
ena. The work presented here extends previous studies 
in two ways: 1. The electrical conductivity and Seeback 
coefficient of the constituent materials are varied inde- 
pendently. It's found that in certain parameter regimes 
the thermoelectric power factor of the composite medium 
is enhanced relative to that of the constituent materi- 
als. 2. The electron and phonon contributions to the 
total thermal conductivity are calculated separately. Al- 
though the constituent materials are assumed to be obey 
the Weidemann- Franz law (W-F) , the composite medium 
does not. Finally, in accordance with previous studies 
[Toj . the figure of merit ZT of the composite medium 
is found to be smaller than that of the high-ZT con- 
stituent material. These conclusions are illustrated with 
numerical results and analytic expressions derived within 
effective medium theory (EMT). 
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where a (r) is the local electrical conductivity, n' e (re' j 
is the electron (phonon) contribution to the total local 
thermal conductivity K (r) (n (r) = k' b (r) + k' (r)) (all 
thermal conductivities evaluated for zero electric field), 
V is the electrostatic potential, and T is the temperature. 
The local Seeback coefficient S (r) is related to a (r) and 
L 12 (r) by: L 12 (r) = S(r)a(r). I assume that a (r) 
and K e (r) obey the W-F law: n e (r) = a (r) Lo^jWhere 
Lo is the Lorenz number. As shown in Ref. (0), the 
effective medium electrical conductivity a, total thermal 
conductivity R, and L 12 satisfy: 
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The brackets indicate an average over disorder configura- 
tions. In addition, the electron and phonon parts of the 
thermal conductivity of the effective medium, R e and k 7 , 
satisfy: 
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Two-component materials are considered here, so that 
local material parameters may assume one of two possible 
values. I further suppose that a\ <C 02, and k 7 i ; 2 "C K e %. 

The figure of merit ZT can be expressed with the On- 
sager number iV (note that N is constrained thermody- 
namically to be less than 1 [n], [l2| ) . 



II. MODEL DESCRIPTION 

The starting point is the linear response description of 
transport for the electrical current j and thermal current 
3Q- 
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In addition to the analysis of EMT, Eq. ([T]) and the 
continuity equations for heat and charge (V • j = 0, 



2 



V • jq =0 [HI) are solved directly for an ensemble of ran- 
domly disordered configurations in 3-d (correlated disor- 
der does not change the results appreciably). The sys- 
tem is discretized into 30 3 sites, and the ensemble size 
is chosen such that the statistical error of the effective 
transport parameters is converged (this typically requires 
about 30 configurations). The error bars on the plots of 
numerical results indicate the statistical uncertainty (one 
standard deviation). 

III. RESULTS 

Power factor - Figs. 1(a), 1(b), and 1(c) show a, L 12 , 
and the power factor S 2 a, respectively, as a function of 
concentration of material 1 (denoted by c). For this 
calculation L{ 2 = L^ 2 . At the percolation threshold 
(c = 2/3), a shows a well-known kink [l4|, while L 12 
is maximized. The origin of this enhancement in L 12 
is discussed in the next section. The enhancement in 
L 12 leads to a peak in the the power factor S 2 a near 
the percolation threshold. The figure shows good agree- 
ment between EMT and numerical results, indicating the 
EMT captures the essential physics of the power factor 
enhancement. Fig. 1(d) shows the Seeback coefficient 
versus concentration when the correct effective medium 
L 12 value is used to calculate S, and when it is assumed 
that L 12 of the composite medium is the same as that of 
the constituent materials. 
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FIG. 1: (a), (b), (c) show the effective conductivity, L 12 , 
and power factor as a function of disorder concentration, (d) 
shows the Seeback coefficient S = P/a as a function of con- 
ductivity (note the log scale of conductivity), using L 12 — 
L 12 , L 12 = L\ 2 . Model parameters are: o\ = 4 X 10 5 
f in3 „ o W i W rl2 rl2 rr~„ 
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To explain the peak in power factor, I solve Eqs. 



at the percolation transition (c = 2/3), and expand the 
solution in the small parameters: 

f s = ^; h = ^. (7) 

0-2 K2 

Keeping only leading order terms leads to: 

L" = Ll 2 (s/2fr: + VV,y 1 + ■■■ (9) 

Eq. ([9]) is valid when L3> 2 is n °t significantly (e.g more 
than 100 times) greater than L\ 2 . Eqs. © and© show 
that both <7 and L 12 diverge as /s" 1 ^ 2 at the percolation 
transition. This implies that the power factor ( (-E 12 ) / <x) 

— 1/2 

also diverges as f s at this point, while the Seeback co- 
efficient (L 12 /a) remains bounded within EMT. The di- 
vergence of transport parameters at the percolation tran- 
sition signals a failure of effective medium theory in the 
limit f s —)-0. However the numerics indicate that EMT 
is sufficient for realistic parameters, and can illustrate the 
key points in this work. 

To develop a simple picture for how inhomogeneity 
leads to an enhancement in L 12 , 1 consider Eq. ([IJ for two 
dissimilar materials placed in series (see Fig. ([!])), where 
<Ti <C <72, ki <C K2, and L\ 2 = L^ 2 - L 12 of the composite 
system is given by ^|av=o- The large thermal resis- 
tance of material 1 implies the temperature drop occurs 
mostly between sites 1 and 2. This drives a large ther- 
moelectric charge current between sites 1 and 2, which 
in turn induces a voltage at site 2. This voltage drives 
an Ohmic charge current between sites 2 and 3. This 
Ohmic current, derived from the inhomogeneous volt- 
age, is the source of enhanced charge current (and there- 
fore enhanced L 12 ) of the composite system. The same 
scenario occurs in 3-dimensional disordered networks of 
materials (under the same parameter region described 
above) , so that when the system is maximally disordered 
at the percolation transition, L 12 shows the maximum 
enhancement. 

Thermal conductivity - Next I consider the electronic 
and phonon contributions to the thermal conductivity of 
the disordered material. Let k 7 i > k i2 ', this describes a 
system in which one material has good electronic thermo- 
electric properties (i.e. a large Onsager number), but a 
detrimentally large k 7 , while the other material has a low 
Onsager number and a small «~. I assume that both 
and are small (of the same order as ^ £l )- Fig. (J3^i) 
shows K e and function of concentration. There 

is good agreement between numerical results and EMT. 
Fig. (J3Jd) shows that near the percolation threshold, the 
electronic contribution to the thermal conductivity is not 
related to the electrical conductivity via the W-F law. At 
the percolation threshold, the expressions for the total 
thermal conductivity k and its partition into electronic 
and phonon parts K e , k 7 are given as (to linear order in 
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FIG. 2: Cartoon indicating the underlying mechanism for the 
enhancement of L 12 for an inhomogeneous systems consisting 
of 2 resistors, (a) shows the predominant paths for charge 
current, (b) shows the voltage and temperature profiles for 
an applied temperature difference. 
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(10) 
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The ratio a/n e is easily expressed in terms of r = fk/fs'- 

(13) 



a 1 
— = L T— 



This indicates that in inhomogeneous materials, inferring 
K e from a may not always be appropriate near the per- 
colation transition. This is because the effective is a 
convolution of the intrinsic material properties and the 
geometry of the system: The parallel conduction of heat 
current through both electron and phonon channels (in 
contrast to the charge current) results in different con- 
ducting paths for jq and j, and differences in the spatial 
structure of temperature and potential fields. As a result, 
the electronic part of j'q in the inhomogeneous system is 
not exactly correlated with j. 

Figure of merit - The electronic and phonon contri- 
butions to ZT can now be assembled. At percolation, 
the Onsager number and ratio of electronic to phonon 
thermal conductivity are: 
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Plugging the above into Eq. ((6]) shows that ZT of the 
composite medium is always smaller than that of the 
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FIG. 3: (a) shows the electron and phonon contributions to 
the total thermal conductivity as a function of concentration 
of material 1 (this concentration is denoted by c) . Solid lines 
are EMT result, symbols are numerical results, (b) shows the 
electrical conductivity and electron contribution to thermal 
conductivity (scaled by the value of material 1 versus concen- 
tration (note log scale). Vertical dashed lines indicate region 
that is plotted in (d). (c) shows the figure of merit ZT as 
a function of concentration. It generally decreases from the 
high material value of ZT — 1. (d) is a zoom- in of (b) near 
percolation (on a linear scale), showing the deviation from the 
W-F law of the composite medium 



high-ZT material constituent. Fig. (j3J;) shows ZT as 
a function of concentration. It is seen that under the fa- 
vorable parameter set considered here, the high ZT value 
is maintained for a substantial amount of low ZT mate- 
rial doping. Generally ZT shows a decrease by percola- 
tion, and rises rapidly to the high value just above this 
threshold. 

I note that this model ignores interface scattering con- 
tributions to <7, k, and L 12 . If these contributions lower k 
and a more drastically than L 12 , it would imply an over- 
all enhancement in ZT 3]. A more complete theory re- 
quires a combination of bulk modeling as presented here, 
combined with more microscopic modeling of materials 
interfaces. Nevertheless the model demonstrates that dis- 
order in the diffusive regime can explain the enhancement 
in power factor, and can alter the relationship between 
electrical and thermal conductivity of electrons. These 
are both important considerations in the interpretation 
of experimental results, and in assessing the viability of 
materials for improved thermoelectric efficiencies. 

I gratefully acknowledge very helpful conversations 
with Fred Sharifi, who introduced me to this problem. 
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